THEP 01/13 
TTPOl-19 
hep-ph/0108197 
August 2001 

The Cross Section of e~^e~ Annihilation into Hadrons of Order a'^n'j 

in Perturbative QCD 



P. A. Baikov 

Institute of Nuclear Physics, Moscow State University, Moscow 119899, Russia 



o 
o 

(N 
< 

m 
(N 



> 

I 

00 

O 



Oh! 
I 



X 



K. G. ChetyrkinQ 

Fakultdt fiir Physik Albert-Ludwigs-Universitdt Freiburg, D-79104 Freiburg, Germany 

J. H. Kiihn 

Institut fiir Theoretische Teilchenphysik, Universitdt Karlsruhe, D-76128 Karlsruhe, Germany 



We present the first genuine QCD five-loop calculation of the vacuum polarization functions: analyt- 
ical terms of order q^Wj to the absorptive parts of vector and scalar correlators. These corrections 
form an important gauge-invariant subset of the full 0{ai) correction to e^e~ annihilation into 
hadrons and the Higgs decay rate into hadrons respectively. They discriminate between different 
widely used estimates of the full result. 
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The value of the strong coupling as is one of the fun- 
damental constants of Nature, comparable in its impor- 
tance to the electromagnetic fine structure constant 
or the weak coupling aweak, characterizing the coupling 
strength of the charged W boson to lefthanded fermions. 
A precise measurement of as is important for the predic- 
tion of numerous observables in the framework of Quan- 
tum Chromodynamics (QCD). The comparison between 
the values determined with the methods of perturbative 
QCD (pQCD) at high energies and from lattice simula- 
tions respectively gives insight in to the quality of approx- 
imation methods and the validity of the overall frame- 
work. The combination of precise values of the three 
couplings as, aom and ctweak in the context of Grand 
Unified Theories leads, last but not least, to important 
information on the structure and the particle content of 
these theories beyond the Standard Model. 

From the theoretical viewpoint the cross section 
for electron-positron annihilation into hadrons (or its 
variants like the .Z-boson decay rate into hadrons or 
the semileptonic branching ratio of the r-lepton into 
hadrons) can be reliably predicted in pQCD and thus 
leads to one of the "gold-plated" as determinations. 

Nontrivial QCD effects in these observables start at or- 
der a^ and their early evaluation has laid the ground 
for the subsequent developments, both on the experimen- 
tal and theoretical side (for a review see, e.g. j|]). Dur- 
ing the past years, in particular through the analysis of 
Z decays at LEP and of t decays, an enormous reduction 



of the experimental uncertainty (down to 3%) has been 
achieved with the perspective of a further reduction by a 
factor of four at a future linear collider |^ . 

Inclusion of the 0{af) corrections 0| is thus manda- 
tory already now. Quark mass effects as well as correc- 
tions specific to the axial current Q must be included 
for the case of ^'-decays. The remaining theoretical un- 
certainty from uncalculated higher orders is at present 
comparable to the experimental one Q. This has led to 
numerous estimates of the C(Q:f ) terms based on princi- 
ples like that of "minimal sensitivity" (PMS) , or of fastest 
apparent convergence (FAC) or by invoking the large (3o 
limit ("naive nonabelization" (NNA)) etc. ||^, ^, ^. In 
fact, already now the extraction of as from r-decays de- 
pends crucially on the corresponding estimates (see, e.g. 

The new calculational methods used in |pj, in partic- 
ular recurrence relations based on the method of inte- 
gration by parts |10 , had opened the door also to the 
four-loop calculation. Nevertheless, it took more than 
ten years of refinement of theoretical tools and develop- 
ment of the algebraic programs to move from the three 
to the four-loop calculation Q. 

It should be noted that much of this effort went into de- 
velopment of the algebraic program MINCER writ- 
ten in FORM , which implements the integration- by- 
parts algorithm to evaluate massless three loop propaga- 
tors. An important part of the calculation, the gener- 
ation of diagrams and their ultraviolet renormalization. 
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was performed manually. In this context it is essential 
that one initially encounters the absorptive parts of four- 
loop amplitudes whose evaluation is reduced to that of 
three- loop propagators, employing the R* operation [ p^ . 
That latter step was also performed manually — an error- 
prone and time-consuming procedure difficult to imple- 
ment in computer algebra, hence presenting a real obsta- 
cle for the calculation. 

Recently the order 0{af) result has been rederived in 
a general covariant gauge in a completely automatic way 
iQ. This latter work was based on a general formula 
which explicitly expresses the result in terms of well- 
known renormalization constants and three-loop mass- 
less propagators, thus resolving the complicated combi- 
natorics of ultraviolet and infrared subtractions. This 
formula is valid for any number of loops and is, therefore, 
an important ingredient of the present 0{a^) calculation. 

The calculational effort for a full evaluation of R{s) in 
0{a'^) is enormous and exceeds the presently available 
computer resources by perhaps two orders of magnitude. 
It is for this reason that we limit the present calculation 
to a gauge invariant subset, namely the terms of order 
a^nj, where n/ denotes the number of fermion flavours. 
The terms of order afn^ (and, in fact, all terms of order 
as(asn/)") have been obtained earlier by summing the 
renormalon chains [ p5[ . These give interesting insight 
into the behavior of perturbation series in high orders and 
have been confirmed by the present work. However, they 
are numerically small and not directly sensitive to the 
nonabelian character of QCD, in contrast to the terms of 
order at 



The motivation of the present work is threefold. First: 
presenting the first 0{a'^) results in pQCD demonstrates 
that the methods proposed in 16, |l^ are indeed 
suited to obtain genuine QCD results in five loop ap- 
proximation. Second: the results per se are an impor- 
tant ingredient of R{s). And third: when comparing the 
exact a'^rij results with the estimates based on various 
optimization procedures one obtains important insights 
into the quality of different approaches, confirming some 
and refuting others. 

Below we will present the results for the vector correla- 
tor and the corresponding quantity R{s) discussed above, 
and for the semileptonic branching ratio of the r lepton. 
The same methods are applied for the scalar correlator. 
In this case the results are of relevance for the Higgs de- 
cays into quarks, for the QCD sum rule determinations of 
the strange quark mass and, last not least, they provide 
a second independent check of the "optimization proce- 
dures" . 

Let us now briefly describe the methods used for the 
practical calculations. Using the techniques from |Q and 
the known renormalization constants up to order the 
four loop propagator type massless integrals remain to 
be evaluated. They were calculated following the method 
described in p6|, p^. It is based on the integration by 



part (IBP) technique however with important mod- 
ifications. The integral is to be reduced to a sum of 
irreducible ( "master" ) integrals with coefficients that are 
known to be rational functions of the space-time dimen- 
sion D. However, in contrast to the original approach, 
where these coefficients were calculated by a recursive 
procedure, they are in the current approach obtained 
from an auxiliary integral representation]!^ . 

For three-loop vacuum integrals these representations 
can be solved explicitly in terms of Pochhammers sym- 
bols . The structure of four- loop propagators is much 
more involved and an explicit solution is not available. In 
the present case the following method is used: for mass- 
less propagators the coefficients of the master integrals 
should be rational functions of space-time dimension D 
only. Their integral representation is expanded in the 
limit 1/D — > 0. Calculating sufficiently many terms in 
this expansion, the original ZJ-dependence can be recon- 
structed. This approach is essentially different from 1/D 
expansion for the very Feynman integrals, recently pro- 
posed in Ref. ||l^ . In our case only the coefficients of the 
master integrals (rational functions in D!) are expanded 
and the original integral can be reconstructed exactly. 

In practice we proceed as follows: First, the set of 
irreducible integrals involved in the problem was con- 
structed, using the criterion irreducibility of Feynman 
integrals jl^. Second, the coefficients multiplying these 
master integrals were calculated in the 1/D — ^ expan- 
sion. This part was performed using the parallel version 
of FORM ||l^ running on an 8-alpha-processor-SMP- 
machine with disk space of 350 GB. The calculations took 
approximately 500 hours in total. (The time includes 
~ 300 hours spent on the single diagram presented on 
Fig. la, which demonstrates that trivial parallelization 
by assigning different diagrams to different machines is 
not feasible.) Third, the exact answer was reconstructed 
from results of these expansions. Extensive tests were 
performed. Dealing with four-loop integrals, the result 
cannot have the singularities stronger then 1/{D — 4)^ 
in D ^ A limit. The contributions corresponding to in- 
dividual master integrals have poles up to 1/(13 — 4)^ 
which cancel in the final expression. Furthermore, sin- 
gular parts of all diagrams were checked by MINCER. 
Previously obtained four-loop results were successfully 
reproduced. 





(b) 



FIG. 1: (a) Diagram whose calculation took ~ 300 hours, (b) 
Diagram contributing to the QED /3-function in 0{a^^n'f). 
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It is convenient to define "Adler functions" as 



(S + Q2)2 



Here, R{s) = 1 + Qs/tt 



[s + Q-'Y ■ 

is tlie famous i?-ratio in 



the massless limit (we suppress the trivial factor J2fQ} 
throughout) and R{s) = 1 + Os/tt + . . . is properly nor- 
malized imaginary part of the scalar polarization func- 
tion. Our results for D,D,R and , R are presented below. 
Let us define (a^ = as/i:) 

oo oc 
oc oo 



where we have set the normalization scale fi^ — or 



to fl 



s for the euclidian and minkowskian functions 



respectively. The results for generic values of /i can be 
easily recovered with the standard RG techniques. Since 
the functions D and D are known to order from 
Refs. g |0) we cite below only the results 
(dots stand for still unknown terms of order a^nf and 
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In (^ and (^ we have distinguished between contribu- 
tions proportional to different color factors, while for (^) 
and (||) their SU(3) values Cp — 4/3, Ca — 3 have been 
specified. In (^ and ^ we differentiate between ^4-terms 
and TT^-ones generated by the analytic continuation from 
spacelike to timelike regions. At last, numerically R and 
R read 

R(s) = l + as + al (1.986 - 0.1153 n/) 
+ al (-6.6369 - 1.2001 n/ - 0.00518 n]) (6) 
+ (0.02152 n) - 0.7974 n) + nj r(4,i) + r(4,o)) , 



R{s) 1 + 5.6667 as + al (35.940 - 1.3586 ri/) 
+ (164.1392- 25.7712 n/-f 0.2590 TiJ) (7) 
+ al (-0.02046 n) + 9.6848 + nj f(4^i) + f(4_o)) . 

The result is easily applied to the semileptonic t decay 
rate. Numerically one finds for the QCD corrections to 
the parton result in the massless limit 
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where we have set fi — and = 3 everywhere except 
for the C'(a^) term. 

Another application of our results is the QED f3- 
function. At four-loops it is known from |2l|. At five- 
loops the only available information is the leading at large 



Uf 



term of order a^ni 



flavours we find a partial 



For QED with Uf electron 
result for the term of order 



^'^n'j), viz. 
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where A = a/(47r),Q; = e^/(47r) and e is the electron 
electric charge. 



4 



TABLE I: Comparison of the results obtained in this paper 
with estimates based of FAC, PMS, NNA and two types of 
Pade approximation methods. 



r-(4,2) 


n-i,2) 




-0.797 


9.685 


exact 


-1.06 
0.74 


7.71 
25.2 
1.00 


|2|](PMS), H (PAM) 
g (PMS, FAC) 
|i (APAM) 
Igll (NNA) 



In Table I our results are compared with predictions 
obtained with the help of various optimization proce- 
dures. To get these numbers we have taken the pre- 
dictions of H, H, H, m, 1^ (where available) for 
n/ — 3,4,5, have subtracted the ©(afn^) contribution 
and fitted the result against the quadratic polynomial in 
n/. For both correlators the predictions based on the 
Principle of Minimal Sensitivity (PMS) or on fastest ap- 
parent convergence (FAC) or, at last, on the Pade Ap- 
proximation Method (PAM) are in quite good agreement 
to the exact results. On the other hand, Asymptotic 
Pade-Approximant Method (APAM) as well as Naive 
NonAbelization (NNA) fail to predict even the rough size 
of the a^nj coefficients. All in all one may conclude that 
there is a good chance that 

• the often used prediction by Kataev and Starshenko 
for the complete 0{a^) contribution to the vector 
correlator will be close to reality; 

• the relatively worser apparent convergence of the 
perturbative series for the scalar correlator as re- 
vealed by the calculation of Rcf. will persist in 
higher orders in agreement with FAC/PMS estima- 
tions of Rcf. ||. 

Finally, we note that the present calculation does not 
include all possible topological types of diagrams appear- 
ing in order a^. Nevertheless, our experience shows that 
these new topologies can be solved in the same way, given 
sufficient computer resources. Work in this direction is 
in progress. 
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